% calcs in NCFMASCrO. Starting guesses provided at the bottom of the file
% Using the model of Jennings & Holland (2015) Journal of Petrology



axfile NCFMASCrO   % Name of AX file to use

infolevel 1

% ignore  liq    % for phases to ignore

printguessform yes
printbulkinfo yes

incax no

setexcess    % for setting excess phases

calctatp ask 

setdefTwindow yes  800 2000    % default T window
setdefPwindow yes  0.001 45    % default P window

project no      % no projections



% starting guesses, pasted in from log file:

% --------------------------------------------------------
% at P = 21.4, T = 1381, for: liq cpx ol sp g opx  with liq = 0, g = 0
% --------------------------------------------------------
ptguess 21.4 1381
% --------------------------------------------------------
xyzguess ct(L)         0.191421
xyzguess jd(L)         0.226294
xyzguess di(L)         0.181008
xyzguess fa(L)         0.105508
xyzguess fo(L)         0.184662
xyzguess hm(L)       0.00885582
xyzguess ek(L)       0.00325739
xyzguess x(cpx)        0.102183
xyzguess y(cpx)        0.102688  range  0.000 2.000
xyzguess o(cpx)        0.364524
xyzguess n(cpx)        0.106650
xyzguess Q(cpx)      -0.0260445  range -1.000 1.000
xyzguess f(cpx)       0.0118534
xyzguess cr(cpx)      0.0209376
xyzguess x(ol)         0.101181
xyzguess x(sp)         0.199508
xyzguess f(sp)        0.0152722
xyzguess cr(sp)       0.0755746
xyzguess x(g)          0.131305
xyzguess c(g)          0.119368
xyzguess f(g)        0.00618964
xyzguess cr(g)        0.0731156
xyzguess x(opx)       0.0969867
xyzguess y(opx)        0.205071  range  0.000 2.000
xyzguess c(opx)       0.0744721
xyzguess Q(opx)      -0.0465413  range -1.000 1.000
xyzguess f(opx)       0.0274898
xyzguess cr(opx)      0.0160842
% --------------------------------------------------------



seta no      % set activity
setiso no    % set isopleths

pseudosection yes

% Bulk composition (converted to mol%):
% ------------------------------------------------
%                SiO2  Al2O3   CaO    MgO   FeO   Na2O     O  Cr2O3  
  setbulk yes  38.494  1.776 2.824 50.566 5.886  0.250  0.096 0.109    % KLB-1 Davis et al 2009, 0.3wt% Fe2O3
% ------------------------------------------------

setmodeiso ask 
zeromodeiso ask 


calcg no
dogmin no

calcsdnle no   % CALCULATE UNCERTAINTY

smath no
drawpd yes

modebox no

xinsmath no
tabsmath no


% script to give melt composition in wt%:
%         SiO2 Al2O3   CaO    MgO   FeO   Na2O      O   Cr2O3
convliq      1     0     0      0     0      0      0    0     % SiO2
convliq      0     1     0      0     0      0      0    0     % Al2O3
convliq      0     0     1      0     0      0      0    0     % CaO
convliq      0     0     0      1     0      0      0    0     % MgO
convliq      0     0     0      0     1      0      0    0     % FeO
convliq      0     0     0      0     0      1      0    0     % Na2O
convliq      0     0     0      0     2      0      1    0     % Fe2O3
convliq      0     0     0      0     0      0      0    1     % Cr2O3






*

% storage area, useful for various starting guesses

% --------------------------------------------------------
% at P = 21.4, T = 1381, for: liq cpx ol sp g opx  with liq = 0, g = 0
% --------------------------------------------------------
ptguess 21.4 1381
% --------------------------------------------------------
xyzguess ct(L)         0.191421
xyzguess jd(L)         0.226294
xyzguess di(L)         0.181008
xyzguess fa(L)         0.105508
xyzguess fo(L)         0.184662
xyzguess hm(L)       0.00885582
xyzguess ek(L)       0.00325739
xyzguess x(cpx)        0.102183
xyzguess y(cpx)        0.102688  range  0.000 2.000
xyzguess o(cpx)        0.364524
xyzguess n(cpx)        0.106650
xyzguess Q(cpx)      -0.0260445  range -1.000 1.000
xyzguess f(cpx)       0.0118534
xyzguess cr(cpx)      0.0209376
xyzguess x(ol)         0.101181
xyzguess x(sp)         0.199508
xyzguess f(sp)        0.0152722
xyzguess cr(sp)       0.0755746
xyzguess x(g)          0.131305
xyzguess c(g)          0.119368
xyzguess f(g)        0.00618964
xyzguess cr(g)        0.0731156
xyzguess x(opx)       0.0969867
xyzguess y(opx)        0.205071  range  0.000 2.000
xyzguess c(opx)       0.0744721
xyzguess Q(opx)      -0.0465413  range -1.000 1.000
xyzguess f(opx)       0.0274898
xyzguess cr(opx)      0.0160842
% --------------------------------------------------------

% --------------------------------------------------------
% at P = 2.1, T = 1145, for: liq pl cpx ol cm opx  with liq = 0, cm = 0
% --------------------------------------------------------
ptguess 2.1 1145
% --------------------------------------------------------
xyzguess ct(L)         0.163583
xyzguess jd(L)         0.145183
xyzguess di(L)         0.165338
xyzguess fa(L)        0.0920349
xyzguess fo(L)         0.153798
xyzguess hm(L)        0.0184116
xyzguess ek(L)       0.00638987
xyzguess Ca(pl)        0.723633
xyzguess x(cpx)        0.102123
xyzguess y(cpx)       0.0483697  range  0.000 2.000
xyzguess o(cpx)        0.184815
xyzguess n(cpx)       0.0403316
xyzguess Q(cpx)      -0.0259273  range -1.000 1.000
xyzguess f(cpx)       0.0111457
xyzguess cr(cpx)      0.0390744
xyzguess x(ol)         0.101273
xyzguess x(cm)         0.427073
xyzguess f(cm)        0.0237435
xyzguess cr(cm)        0.830201
xyzguess x(opx)        0.100628
xyzguess y(opx)        0.115361  range  0.000 2.000
xyzguess c(opx)       0.0843996
xyzguess Q(opx)      -0.0463361  range -1.000 1.000
xyzguess f(opx)       0.0384826
xyzguess cr(opx)      0.0337310
% --------------------------------------------------------

% --------------------------------------------------------
% at P = 0.0, T = 1116, for: liq pl cpx ol sp opx  with liq = 0
% --------------------------------------------------------
ptguess 0.0 1116
% --------------------------------------------------------
xyzguess ct(L)         0.138751
xyzguess jd(L)         0.123225
xyzguess di(L)         0.152718
xyzguess fa(L)        0.0894335
xyzguess fo(L)         0.161687
xyzguess hm(L)        0.0216911
xyzguess ek(L)       0.00616605
xyzguess Ca(pl)        0.732671
xyzguess x(cpx)        0.101558
xyzguess y(cpx)       0.0360303  range  0.000 2.000
xyzguess o(cpx)        0.174030
xyzguess n(cpx)       0.0352314
xyzguess Q(cpx)      -0.0256647  range -1.000 1.000
xyzguess f(cpx)       0.0108576
xyzguess cr(cpx)      0.0342950
xyzguess x(ol)         0.101160
xyzguess x(sp)         0.447936
xyzguess f(sp)        0.0248720
xyzguess cr(sp)        0.872106
xyzguess x(opx)        0.100585
xyzguess y(opx)       0.0980281  range  0.000 2.000
xyzguess c(opx)       0.0855804
xyzguess Q(opx)      -0.0460593  range -1.000 1.000
xyzguess f(opx)       0.0392659
xyzguess cr(opx)      0.0295954
% --------------------------------------------------------

% --------------------------------------------------------
% at P = 41.9, T = 1640, for: liq cpx ol g opx  with liq = 0, opx = 0
% --------------------------------------------------------
ptguess 41.9 1640
% --------------------------------------------------------
xyzguess ct(L)         0.126479
xyzguess jd(L)        0.0621998
xyzguess di(L)         0.232302
xyzguess fa(L)         0.144751
xyzguess fo(L)         0.318568
xyzguess hm(L)        0.0196377
xyzguess ek(L)       0.00325383
xyzguess x(cpx)       0.0870305
xyzguess y(cpx)       0.0629732  range  0.000 2.000
xyzguess o(cpx)        0.569641
xyzguess n(cpx)       0.0746353
xyzguess Q(cpx)      -0.0201195  range -1.000 1.000
xyzguess f(cpx)       0.0178205
xyzguess cr(cpx)     0.00686655
xyzguess x(ol)         0.102464
xyzguess x(g)          0.124677
xyzguess c(g)          0.112268
xyzguess f(g)         0.0277346
xyzguess cr(g)        0.0656923
xyzguess x(opx)       0.0826396
xyzguess y(opx)        0.141426  range  0.000 2.000
xyzguess c(opx)       0.0999616
xyzguess Q(opx)      -0.0331937  range -1.000 1.000
xyzguess f(opx)       0.0408887
xyzguess cr(opx)     0.00630170
% --------------------------------------------------------

